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1.  Introduction 


This  work  is  part  of  a  cooperative  research,  development,  test  and  evaluation  program  that  brings  scien¬ 
tists  and  engineers  from  various  English-speaking  countries  together  for  collaborative  studies  to  improve 
technology  to  solve  technical  problems.  The  Technical  Cooperation  Program  (TTCP)  identified  Computa¬ 
tional  Fluid  Dynamics  (CFD)  as  a  critical  technology.  The  current  TTCP  work  project  WTP-2  KTA  2-12, 
“Application  of  CFD  to  the  prediction  of  Missile  Body  Vortices”  extends  on  knowledge  gained  by  KTA 
2-9  and  focuses  on  the  ability  of  the  three-dimensional  Navier-Stokes  equations  to  predict  fiowfields  about 
high  length- to- diameter  bodies  at  moderate  angles  of  attack  (8®  <  o'  <  14®)  for  transonic  and  supersonic 
Mach  numbers. 

This  report  documents  this  authors’s  computational  results  for  the  five  test  cases.  The  five  test  Ccises 
are  as  follows. 

Case  1:  Mach  1.45  at  cr  =  14®,  Red  =  616,666. 

Case  2:  Mach  1.8  at  o*  =  14®,  Rej^  =  616,666. 

Case  3:  Mach  2.5  at  cr  =  14®,  Re^  =  1,233,333. 

Case  4:  Mach  3.5  at  cr  =  8®,  Rcd  =  1,233,333. 
and  Case  5:  Mach  3.5  at  cr  =  14®,  Ren  =  1,233,333. 

The  experimental  body  is  13  diameters  long  with  the  diameter  of  the  cylindrical  afterbody  of  3.7  inches. 
Laminar  and  turbulent  computations  are  shown  with  comparisons  to  experimental  data.  The  experimental 
validation  data  were  provided  by  the  Defence  Research  Agency  (UK).  The  computations  were  performed 
with  an  existing  Navier-Stokes  code  FDL3DI  developed  in  the  Wright  Laboratory.  The  code  features  two 
kinds  of  spatial  discretization  which  were  used  for  the  present  studies,  one  having  the  central  differences 
for  both  inviscid  and  viscous  terms  and  the  other  having  an  upwind-biased  discretization  for  the  inviscid 
terms  and  the  viscous  terms  are  central  differenced.  The  code  includes  the  k-e  Turbulence  model  which  was 
used  for  simulating  turbulence  effects  in  the  flowfield.  Grid  resolution  study  was  conducted  to  determine 
important  areas  in  flow  regimes  where  high  grid  densities  are  important  to  resolve  the  vortices  and  shocks 
and  also  for  accurately  predicting  the  surface  pressure  and  pitot  pressure  in  the  flow  field.  In  addition  to 
the  computation  of  five  test  cases  on  grids  that  were  considered  to  be  optimum  for  each  of  the  cases,  a 
standard  grid  was  generated  by  members  of  the  Wright  Laboratory  to  be  used  by  all  the  participants  for 
comparing  results  with  different  CFD  codes  and  turbulent  models.  The  standard  grid  wais  computed  only 
for  Case  3  in  the  present  study. 

At  supersonic  Mach  numbers  the  fiowfields  around  missile  body  configurations  at  angle- of-attack  have 
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boundary  layer  separation  leading  to  vortical  structures.  The  crossflow  separation  on  the  leeward  side 
is  characterized  by  the  crossflow  separation  line  and  the  strength  of  the  vortex  which  develops  from  the 
separation.  Prediction  of  these  quantities  by  comput atonal  methods  is  highly  dependent  on  the  available 
turbulent  models  and  their  improvements. 

There  have  been  a  number  of  numerical  studies  in  recent  years  dealing  with  the  problem  of  vortical  flows 
for  ogive-cylinder  configurations  at  angle  of  attack  [1,  2,  3,  4].  The  importance  of  radial  grid  resolution  for 
the  attached  turbulent  viscous  layers  and  the  radial  and  circumferential  resolution  of  the  leeward  vortices 
was  noted  by  Degani,  et  al.  [1].  Their  work  used  the  algebraic  eddy-viscosity  model  and  modifications 
were  made  to  the  turbulence  model  to  properly  evaluate  the  viscous-layer  scale  length  under  the  leeward 
vortex  structure.  The  computational  work  of  Borrel  et  al.  [2]  also  used  the  algebraic  eddy-viscosity  model 
and  noted  the  problem  of  determining  the  mixing  length  where  there  were  pronounced  turbulent  effects. 
The  work  of  Hsieh,  et  al.  [3]  used  the  algebraic  turbulence  model  modified  by  Degani,  et  al.  [1]  and  found 
that  the  discrepancy  between  calculation  and  experiment  increased  at  10  and  15  degrees  incidences  for  a 
three-inch  diameter  body  at  freestream  Mach  number  of  3.5.  The  work  of  Moran,  et  al.  [4]  also  used  the 
algebraic  turbulent  model  and  made  modifications  to  account  for  compressibility  and  pressure  gradients 
and  also  altered  the  two  coefficients  in  the  outer  region  to  better  match  the  law  of  the  wake.  Computations 
were  performed  as  a  participant  of  the  KTA  2-12  by  Housh  [5]  using  a  TVD  scheme  and  the  Baldwin-Barth 
turbulence  model. 

In  view  of  the  importance  of  turbulence  models  in  modeling  these  flows,  the  present  work  used  the 
two-equation  k  —  e  Turbulence  model  and  compared  it  to  the  laminar  solution  results.  The  study  contains 
numerical  simulations  of  the  Navier-Stokes  equation  with  and  without  turbulent  models  for  supersonic  flows 
at  moderate  angles  of  attack  for  which  detailed  experimental  data  is  available.  Two  different  numerical 
spatial  discretizations  were  used  to  compare  solutions  to  determine  if  there  were  any  issues  pertaining  to 
algorithm  development. 

The  objectives  of  the  current  KTA  2-12  are  to  assess  the  current  capabilities  of  Navier-Stokes  solvers 
for  the  prediction  of  weapon  flowfields  and  to  promote  technical  interchange  among  CFD  practitioners 
and  experimental  personnel  in  the  areas  of  grid  generation,  solution  algorithms,  turbulence  modeling,  code 
performance,  validation,  and  data  visualization.  The  objectives  will  be  accomplished  by  comparing  the 
results  obtained  from  Navier-Stokes  computations  against  experimental  data,  evaluating  the  accuracy  of 
the  predictive  technology  and  identifying  priorities  for  future  development. 
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2.  Analysis 


This  section  gives  details  of  fluid  dynamic  equations  and  turbulence  model  of  an  existing  computer  code 
FDL3DI  developed  in  the  Wright  Laboratory.  The  governing  equations  were  taken  to  be  the  unsteady  three- 
dimensional  compressible  mass-averaged  turbulent  Navier-Stokes  equations,  written  in  nondimensional 
variables  and  expressed  notationally  in  the  following  strong- conservation  form 


dQ  d 


Here  t  is  the  time,  i,  rj,  C  the  surface  fitted  computational  coordinates,  Q  the  vector  of  dependent  variables, 
F,G,  H  the  inviscid  flux  vectors,  Fv,Gv,  Hv  the  viscous  flux  vectors,  and  S  the  vector  source  term.  Effects 
of  turbulence  were  accounted  for  by  specifying  a  turbulent  Prandtl  number  Pr«  =  0.90,  and  by  incorporating 
a  two-equation  model  for  the  turbulence  kinetic  energy  k  and  the  turbulence  dissipation  rate  e.  The  k-e 
turbulence  model  is  that  of  Jones  and  Launder  [6,  7]  for  which  the  generalized  formulation  was  written 
by  Gerolymos  [8].  The  model  incorporates  low  Reynolds  number  terms  to  account  for  near- wall  effects. 
The  turbulence  model  is  implemented  with  compressibility  correction  of  Sarkar,  et  al.  [9]  for  more  accurate 
representation  of  supersonic  flows. 

With  this  formulation,  the  vector  of  dependent  variables  is  given  as 

P 

pu 
pv 

pw  (2.2) 

pE 
pk 
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and  the  vector  fluxes  by 


-7 


pu 
puU  + 
pvU  +  ^yP 
pwU  +^zP 
pEU  +  ^xiWiP 
pkU 
pcU 

pV 

puV  +  TfiP 
PVV  +  TJyP 

pwV  +  r)zp 
pEV  +  r)xiUip 
pkV 
peV 

pW 

puW  +  Cxp 
pvW  +  CyP 
pwW  +  (izP 
pEW  +  GiUiP 
pkW 


J 


^Xi  7t2 
^Xi'^iS 

ixiiujTij  -  qi) 


4 


0 


Vxi'^i2 
'Oxi'^iZ 

Vx.i‘>^jTij  -  qi) 

0 

Cxi^vi 
CxiTi2 
Cxi'^iZ 

Cxi{'iJ>jTij  — 

ifJ'  +  ^^/^k)Cxi'^f^ 

{fi  +  ^i/(T^)(^s:i 

with  the  vector  source  term 

0 
0 
0 

0  (2.9) 

0 

-[Pk-pc{l  +  M^)-VD] 

^{C^\Pk^/k  ~  Cc2ft2P^^ /k  +  E^) 

where 


u  =  ^t 

(2.10) 

V  =  Tit  +  TJxiUi 

(2.11) 

W  =  Ct  +  QiUi 

(2.12) 

T  1 

1)M2,  +  2  + 

(2.13) 

In  the  preceding  expressions,  u,v,w  zxe  the  Cartesian  velocity  components,  p  the  density,  p  the  pressure, 
and  T  the  temperature.  All  length  scales  have  been  nondimensionalized  by  the  diameter  of  the  cylinderical 
afterbody  D  and  dependent  variables  have  been  normalized  by  freestream  values  except  for  p,  jfe,  and  e 
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which  have  been  nondimensionalized  by  respectively.  Components  of  the  heat  flux 

vector  and  stress  tensor  may  be  expressed  as 
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and  the  turbulence  parameters  are 
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=  exp 


/,2  =  1-0.3  exp(-i2e2) 

C;,  =  0.09,  Cel  =  1.44,  Ce2  =  1.92 
cTfc  =  1.0,  <Je  =  1-3. 

The  Sutherland  law  for  the  molecular  viscosity  coefficient  fi  and  the  perfect  gas  relationship 


(2.15) 


(2.22) 


(2.24) 


(2.26) 


were  also  employed,  and  Stokes’  hypothesis  for  the  bulk  viscosity  coefficient  has  been  invoked. 

The  coefficients  of  the  force  components  given  by  the  experiment  are  the  Cx  and  Cz  which  are  defined 
as  follows; 


C force  — 

The  coefficient  of  moment  is  defined  as  follows: 


Force 
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^  _  Moment 

^moment  —  TT"!  t7o"T7^ 

O.bpooV^AD 

The  moments  are  measured  about  the  nose  tip  (X  =  0). 

The  coefficients  of  lift  and  drag  are  defined  by: 

Cl  =  Czcosa  —  CxSincr 


(2.29) 


(2.30) 


Cd  =  Czsincr  +  Cxcoscr  (2.31) 

In  the  above,  the  reference  length  D  is  the  diameter  of  the  body  which  is  3.7  in  and  the  Area  A  is  given 
by  ^ 

Except  for  the  axial  force  Cx,  all  the  other  comparisons  are  got  by  integrating  the  pressures  and  shear 
stresses  on  the  body.  The  axial  force  experimental  results  are  disturbed  by  the  sting  mounting  and  in  the 
computation  the  freestream  static  pressure  is  applied  over  the  base  area. 

The  components  of  the  stress  tensor  written  in  the  non-dimensional  variables  is  given  by, 
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(2.36) 
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ryz-rzy-—{p  +  pz){—  +  —) 

(2.37) 

The  components  of  force  in  the  transformed  coordinates  are  given  by, 
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The  mass  averaged  Navier-Stokes  equations  were  solved  using  two  kinds  of  spatial  discretization  both 
based  on  the  implicit  Beam  and  Warming  algorithm.  The  central  difference  version  implemented  into  the 
code  is  2nd  order  accurate  spatially  and  uses  a  blended  2nd/4th  order  Jameson-type  damping.  The  other 
kind  of  spatial  discretization  which  used  the  Roe  Approximate  Riemann  solver  for  the  inviscid  terms  which 
has  3rd  order  spatial  accuracy  and  central  differences  for  the  viscous  terms.  The  Navier-Stokes  equations 
solved  by  the  Beam  and  Warming  algorithm  were  used  successfully  in  analyzing  vortex  breakdown  above 
a  pitching  delta  wing  by  Visbal  [10]  and  turbulent  cylinder  juncture  flows  by  Rizzetta  [11].  The  Roe  solver 
was  used  to  study  the  turbulent  interactions  due  to  cylinder/offset  flare  juncture  by  Gaitonde,  et  al.  [12]. 


2.1  Boundary  Conditions 

Since  supersonic  conditions  were  at  both  inflow  and  outflow,  at  the  inflow,  supersonic  conditions  were  ap¬ 
plied  and  at  outflow  extraplolation  condition  was  used.  Along  the  farfield  boundary,  freestream  conditions 
were  given.  For  the  body  surfaces,  no  slip  was  applied,  zero  normal  pressure  gradient  and  an  adiabatic 
wall  was  used.  For  the  turbulent  case,  homogeneous  conditions  for  k  and  e  were  specified. 

A  line  of  singularity  in  the  coordinate  transformation  lies  along  the  longitudinal  axis  and  extends  from 
the  tip  of  the  ogive  to  the  upstream  farfield  boundary.  The  flow  variables  along  this  line  (i=l)  were  simply 
set  equal  to  the  average  of  the  flow  variables  on  the  i=2  surface.  The  flow  is  assumed  to  have  bilateral 
symmetry  with  respect  to  the  x-y  plane;  therefore  w=rO  and  ^  =  0  for  all  the  remaining  variables. 

2.2  Experimental  Data  and  Conditions  of  Numerical  Simula¬ 
tion 

The  experimental  data  were  provided  by  the  Defence  Research  Agency,  UK.  The  geometry  consists  of  a 
length  of  13  diameters  as  shown  in  Figure  1.  The  nose  is  given  by  the  equation: 

r{x)ID  =  -0.002615(a:/D)3  -  0.03986(a:/D)2  +  0.30984(x/D)  (2.41) 

Fig.  2  shows  the  body  and  the  computational  grid  used  for  computing  the  conditions  of  Case  1.  For 
this  case,  the  asymmetricity  in  the  grid  allows  for  better  capture  of  vortices  on  the  leeward  side  and  the 
strong  shock  close  to  the  body  on  the  windward  side.  The  grid  was  generated  using  the  GRIDGEN3D  [13]. 

The  five  test  cases  chosen  for  this  study  are  the  free  stream  conditions  shown  in  the  Table  2.1.  Priorti- 
zation  of  the  cases  as  determined  by  the  TTCP  workgroup  is  shown  in  Table  2.1. 
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Table  2.1:  Free  Stream  Conditions 


Case 

No. 

Priority 

No. 

Mach 

No. 

Angle  of 

Attack,  deg. 

Reo 

Total 

press.,  N/m^ 

Total 

temp.,  K 

1 

5 

1.45 

14 

616,666 

44,912 

301 

B 

2 

1.8 

14 

616,666 

47,276 

304 

B 

1 

2.5 

14 

1,233,333 

141,827 

308 

B 

3.5 

8 

1,233,333 

243,133 

315 

B 

3.5 

14 

1,233,333 

243,133 

315 

Table  2.2:  Available  Experimental  Data 


Case  No. 

Surface  Pressure,  Cp 

Pitot  Pressure 

1 

X/D=0.2-14.5,  <i>=0-m 

X/D=8.5  and  11.5 

2 

X/D=0.2-14.5,  <?i=0-180 

X/D=5.5,  8.5  and  11.5 

3 

X/D=0.2-14.5,  <?i=0-180 

X/D=8.5  and  11.5 

4 

X/D=0.2-14.5,  «i=0-180 

X/D=5.5 

5 

X/D=0.2-14.5,  <6=0-180 

X/D=5.5  and  11.5 

The  experimental  data  consist  of  surface  pressure  and  pitot  pressure.  The  surface  pressure  is  at  32 
axial  locations  and  the  pitot  pressure  is  at  axial  locations  as  specified  in  the  Table  2.2.  Although,  the 
computations  were  performed  for  a  length  of  13  diameters,  the  experimental  body  was  longer  and  the 
pressure  coefficients  were  taken  upto  a  length  of  14.5  diameters. 

The  grids  developed  for  this  study  are  of  the  C-type  grids.  Since  the  study  was  pertaining  to  steady 
state  solutions,  only  a  half  grid  in  the  circumferential  direction  was  used.  Except  for  the  grid  for  Case  1 
and  the  standard  grid,  the  C-grid  was  rotated  about  the  longitudinal  axis  for  the  body  to  contruct  the 
three-dimensional  grid.  For  the  grid  for  Case  1  and  the  standard  grid,  the  GRIDGEN3D  [13]  was  used 
to  rotate  the  grid  about  the  longitudinal  axis  for  an  asymmetry  required  for  better  resolution  of  vortical 
structures  on  the  leeward  side  and  the  shock  wave  close  to  the  body  on  the  windward  side. 

The  feature  of  local  time  stepping  was  used  for  the  central  difference  version  of  the  code  to  speed  up 
convergence  to  steady  state.  The  constant  nondimensional  time  steps  employed  for  these  calculations,  as 
dictated  by  stability  requirements  are  0.01  <  (r  <  0.0001.  The  turbulent  calculation  was  started  from  the 
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Table  2.3:  Details  of  grids  considered  in  this  study 


Case  No. 

Mach  No. 

AoA 

Red 

A„ 

Grid  size 

1 

1.45 

14'’ 

616,666 

0.00001 

173  X  99  X  93 

2 

1.8 

14® 

616,666 

0,00003 

150  X  82  X  75 

3 

2.5 

14® 

1,233,333 

0.00003 

100  X  50  X  51 

3 

2.5 

14® 

1,233,333 

0.00003 

150  X  75  X  76 

3 

2.5 

14® 

1,233,333 

0.00001 

121  X  89  X  89 

4 

3.5 

8® 

1,233,333 

0.00003 

150  X  75  X  76 

5 

3.5 

14® 

1,233,333 

0.00003 

150  X  75  X  76 

An  =  Minimum  Normal  Grid  Spacing  normalized  by  Rn  =  3.7m 
Grid  size  is  given  as:  axial  x  radial  x  circumferential 


converged  laminar  solution. 


2.3  Grid  Study 

A  grid  study  was  conducted  to  determine  the  optimum  grid  required  to  compare  the  computational 
results  to  data.  The  grid  study  results  are  shown  only  for  Case  3.  The  grid  sizes  and  minimum  distances 
for  all  three  grids  for  Case  3  and  for  the  other  cases  is  given  in  Table  2.3.  The  laminar  and  turbulent  pitot 
pressure  along  the  vortex  centerline  for  the  three  grids  for  Case  3  is  shown  at  an  X/D  location  of  11.5  in 
Figures  3  and  4.  It  can  be  seen  that  as  the  points  in  the  radial  and  circumferential  direction  are  increased, 
the  vortex  shape  is  better  captured.  The  surface  pressure,  however,  for  all  three  grids  have  small  variations 
with  each  other  (figures  not  shown).  From  the  comparison  of  laminar  and  turbulent  computations  it  can 
also  be  seen  that  the  vortex  core  is  better  captured  in  the  turbulent  case. 
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3.  Computer  Requirements 


The  computational  resources  used  for  computing  the  test  cases  employing  three  different  grids  used  in 
the  present  study  are  summarized  in  this  section.  The  summary  shown  in  Table  3.1  gives  the  computer 
resources  used  for  computational  results  presented  in  this  report. 

The  Data  Processing  Rate  (DPR)  for  the  laminar  cases  for  both  Central  and  Upwind  difference  versions 
are  close  to  each  other.  The  DPR  for  the  turbulent  cases  are  significantly  higher  due  to  the  inclusion  of 
the  two-equation  k-e  turbulence  model.  Part  of  the  reason  for  the  high  DPR  is  because  the  computer  code 
for  the  turbulence  model  is  not  fully  optimized. 

The  Cray  C-90  was  used  for  computation  of  all  test  cases 
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Table  3.1:  Computer  resources  requirements 


Grid  Size 

Type  of  discretization 

Laminar  or  Turbulent 

DPR 

Memory 

100  X  50  X  51 

Central 

Laminar 

2.3 

10  Mw 

100  X  50  X  51 

Central 

Turbulent 

5.4 

10  Mw 

150  X  75  X  76 

Central 

Laminar 

7.7 

33  Mw 

150  X  75  X  76 

Central 

Turbulent 

18.0 

33  Mw 

121  X  89  X  89 

Central 

Laminar 

8.3 

37  Mw 

121  X  89  X  89 

Central 

Turbulent 

20.2 

37  Mw 

121  X  89  X  89 

Upwind 

Laminar 

11.3 

37  Mw 

121  X  89  X  89 

Upwind 

Turbulent 

24.5 

37  Mw 

DPR  =  Data  Processing  Rate 

Grid  size  is  given  as:  axial  x  radial  x  circumferential 
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4.  Discussion  of  Results 


This  section  has  discussion  of  results  for  computations  of  all  five  test  cases  and  the  computations  on 
the  standard  grid  for  the  Case  3.  The  computations  on  the  standard  grid  use  the  two  types  of  spatial 
discretization  -  (a)  central  differences  for  inviscid  and  viscous  terms  and  (b)  upwind  differences  for  inviscid 
terms. 


4.1  Surface  Pressure  Comparisons  for  Case  3  on  Standard  Grid 

Figures  5-  20  show  the  computed  surface  pressure  comparisons  with  experiments.  The  comparisons  are 
for  both  versions  of  the  code,  namely  central  and  upwind  difference  and  for  both  laminar  and  turbulent 
computations.  The  results  are  shown  for  various  axial  stations  of  X/D  =  2.4,  3.5,  4.5,  5.5,  6.5,  7.5,  8.5,  and 
11.5.  For  all  the  results  shown,  the  surface  pressures  obtained  by  central  difference  and  upwind  difference 
are  very  close  to  each  other.  The  laminar  and  turbulent  computations,  however,  are  showing  differences 
after  <j)  of  90*  for  all  stations,  the  variation  becoming  greater  as  X/D  increases.  The  comparison  of  turbulent 
computations  agree  with  the  experiment  better  than  the  laminar  computations  at  all  X/D  locations.  The 
laminar  computations  show  additional  suction  peaks  beyond  of  60*  for  all  locations  denoting  additional 
vortical  structures  in  the  flow  field.  The  turbulent  computation  underpredicts  the  second  suction  peak  at 
X/D  of  4.5,  5.5,  and  8.5  (shown  in  Figures  10,  11,  and  18)  which  can  be  attributed  to  excess  dissipation  of 
the  k-c  Turbulence  model.  The  surface  pressure  beyond  120*  at  X/D  of  8.5  shown  in  Figure  18  shows  poor 
comparison  with  experimental  values.  The  surface  pressure  comparisons  between  the  central  difference  and 
the  upwind  difference  versions  show  very  little  variations  for  both  laminar  and  turbulent  computations. 


4.2  Pitot  Pressure  Comparisons  for  Case  3  on  Standard  Grid 

The  pitot  pressure  contours  predicted  by  computations  and  the  pressure  extracted  along  the  vortex  cen¬ 
terline  are  compared  to  the  pitot  pressure  contours  and  pressure  extracted  along  the  vortex  centerline  of 
the  experiment  in  Figures  21  to  32.  The  computational  results  shown  are  for  both  laminar  and  turbulent 
for  the  central  and  upwind  difference  versions  on  the  standard  grid.  The  pitot  pressure  of  the  computation 
outside  of  the  vortex  is  same  as  the  experiment.  The  shock  wave  at  an  angle  of  about  45*  of  the  circular 
cross-section  for  the  experimental  pitot  pressure  contours  is  better  captured  for  the  turbulent  case  than  for 
the  laminar,  as  can  be  seen  in  Figures  21  and  22.  The  pitot  pressure  along  the  vortex  centerline  shows  the 
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turbulent  computations  to  be  better  than  laminar  shown  in  Figures  23  and  26;  however,  at  X/D  of  11.5, 
the  vortex  pressure  of  the  laminar  case  matches  the  experimental  value  better,  as  seen  in  Figures  29  and 
32.  For  the  laminar  computation,  additional  vortices  can  be  seen  at  X/D  of  11.5  as  seen  in  Figure  27  and 
30. 


4.3  Pitot  Pressure  Comparisons  for  Cases  1  to  5 

The  pitot  pressure  comparisons  and  pressure  along  vortex  centerline  for  Case  1  are  shown  in  Figures  33 
to  38.  Additional  suction  peaks  are  seen  in  the  turbulent  computations  at  X/D  of  8.5  and  11.5  as  seen 
in  Figures  35  and  38.  At  X/D  of  11.5  the  vortex  shows  excessive  dissipation  due  to  the  effects  of  the 
turbulence  model  as  seen  in  Figures  37  and  38. 

The  pitot  pressure  comparisons  and  pressure  along  vortex  centerline  for  Case  2  are  given  in  Figures  39 
to  46.  There  is  no  well-defined  vortex  in  the  computation  at  X/D  location  of  5.5,  hence  no  comparison 
is  made  for  pressure  along  vortex  centerline  with  experiment.  At  X/D  of  8.5,  the  laminar  computation 
compares  with  experiment  better  than  the  turbulent  case  as  can  be  clearly  seen  from  the  comparison  of 
pressures  along  the  vortex  centerline.  At  X/D  of  11.5,  both  laminar  and  turbulent  pitot  pressures  show  the 
vortex  to  be  more  dissipated  than  the  vortex  in  experiment.  However,  the  shape  and  size  of  the  laminar 
vortex  compares  well  with  the  experiment.  The  pitot  pressure  comparisons  for  Case  3  are  shown  in  Figures 
47  to  51.  There  is  no  clear  vortex  center  in  the  computation  at  X/D  of  5.5  hence,  no  pitot  pressure  along 
vortex  centerline  is  shown  in  comparison  to  experiment.  At  X/D  of  11.5,  the  comparison  between  the 
laminar  and  turbulent  computations  and  experiment  are  reasonable  with  the  vortex-core  of  the  laminar 
computation  matching  the  experiment  better  than  the  turbulent  case.  The  pitot  pressure  comparisons  and 
the  pressure  along  the  vortex  centerline  for  Case  4  are  shown  in  Figures  52  to  54.  The  vortex  position 
shown  in  Figures  52,  53  and  54  is  better  predicted  for  the  turbulent  case  than  for  the  laminar  case; 
however,  the  turbulent  model  causes  excess  dissipation.  The  pitot  pressure  comparisons  and  the  pressure 
along  vortex  centerline  for  Case  5  are  shown  in  Figures  55  to  59.  The  comparison  between  the  experiment 
and  computation  at  X/D  of  5.5  is  reasonable  as  shown  in  Figures  55  and  56.  At  X/D  of  11.5,  referring  to 
Figures  57  and  58  the  pitot  pressure  contours  of  the  computations  do  not  show  the  details  of  the  vortex 
structure  as  seen  in  the  experiment.  The  vortex  centerline  pressure,  however,  shows  reasonable  comparison 
with  experiment. 
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4.4  Surface  Pressure  Comparisons  for  Cases  1  to  5 


The  surface  pressure  comparisons  for  Case  1  are  shown  in  Figures  60  to  67.  The  surface  pressure  using 
the  turbulence  model  shows  better  comparison  than  laminar.  There  are  as  many  as  four  suction  peaks 
showing  in  some  laminar  profiles.  The  profile  using  the  turbulent  model  at  X/D  of  4.5,  5.5  and  6.5  shows 
the  excess  diffusivity  of  the  turbulence  model.  The  surface  pressure  comparisons  for  Case  2  are  shown  in 
Figures  68  to  75.  The  surface  pressure  using  the  tubulence  model  compares  to  the  experiment  very  well  at 
all  the  X/D  locations  shown.  The  suction  peak  has  shifted  compared  to  the  experiment  for  the  turbulence 
profile  in  the  case  of  X/D  of  6.5.  The  excess  diffusivity  of  the  turbulence  model  can  be  seen  at  X/D  of  5.5, 
6.5  and  7.5.  The  surface  pressure  comparisons  for  Case  3  are  shown  in  Figures  76  to  83.  The  profiles 
using  the  turbulence  model  shows  better  comparison  with  experiment  than  the  laminar.  The  most  notable 
discrepancy  of  the  turbulent  profile  with  the  experiment  is  at  the  X/D  location  of  7.5;  after  the  suction 
peak  at  ^  =  90®,  the  turbulent  computation  matches  the  laminar  better  than  the  experiment. 

The  surface  pressure  comparisons  for  Case  4  are  shown  in  Figures  84  to  91.  The  turbulent  computation 
profiles  match  the  experiment  very  well.  At  X/D  locations  of  7.5  and  8.5,  there  is,  however,  excess 
dissipation  due  to  the  turbulence  model  which  gives  a  higher  suction  pressure  at  0  =  150®.  The  surface 
pressure  comparisons  for  Case  5  are  shown  in  Figures  92  to  99.  The  turbulent  profiles  compare  with  the 
experiment  better  than  the  laminar  ones.  The  turbulent  profile  at  X/D  of  6.5  mathces  the  laminar  profile 
better  than  the  experiment. 


4.5  Convergence  History  for  Case  3  on  Standard  Grid 

The  convergence  history  of  the  solution  using  the  upwind  version  of  the  code  to  reach  steady  state  is  shown 
in  Figures  100  and  101.  For  the  laminar  case  the  L2norm  is  reduced  to  six  orders  of  magnitude  and  for 
the  turbulent  case  convergence  is  assumed  when  the  normalized  surface  pressure  (averaged  over  the  entire 
surface  of  the  body)  is  shown  to  reach  steady  state. 

4.6  Forces  and  Moments  for  Case  3  on  Standard  Grid 

The  axial  and  radial  components  of  the  force  and  the  moment  about  y  was  computed  for  Case  3  using  both 
the  central  and  upwind  difference  versions  of  the  code.  The  Table  4.1  gives  the  comparison  of  forces  and 
moments  for  Case  3  on  the  standard  grid. 

Referring  to  Table  4.1,  it  can  be  seen  that  the  predictions  by  turbulent  computations  made  with  both 
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Table  4.1:  Comparison  of  Forces  and  Moments  with  Experiment  for  Central  and  Upwind  Difference  Versions 


Refer  Figure  1  for  sign  convention 

Experimental  values  are  interpolated  from  data  at  a  of  12.16®  and  14.20° 

the  central  and  upwind  difference  versions  match  the  experimental  Cx,  Cz,  Cl,  Cd,  and  Cm  better  than 
the  predictions  by  laminar  computations. 
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5.  Summary  and  Conclusions 


As  part  of  TTCP  work  project,  the  KTA  2-12  “Application  of  CFD  to  the  prediction  of  missile  body 
vortices,”  computations  were  performed  to  assess  the  predictive  capability  of  full  Navier-Stokes  equations 
with  turbulence  model  for  two  different  types  of  discretization  for  the  convective  terms.  An  existing 
Navier-Stokes  code  FDL3DI  developed  in  the  Wright  Laboratory  was  used  to  perform  the  computations. 
In  the  code  the  mass  averaged  Navier-Stokes  equations  were  solved  with  the  implicit  Beam  and  Warming 
algorithm  using  (a)  the  central  difference  for  both  convective  and  viscous  terms  with  the  required  numerical 
damping  and  (b)  the  upwind  difference  discretization  for  convective  terms  and  central  differences  for  viscous 
terms  without  explicit  numerical  damping.  The  turbulence  model  which  is  part  of  FDL3DI,  chosen  for  this 
study  was  the  k-e  Turbulence  model.  The  turbulence  model  included  compressiblity  correction  and  the 
low  Reynolds  number  terms  to  include  near-wall  effects.  All  computational  results  were  compared  to  the 
experimental  data  provided  by  the  Defence  Research  Agency  (DRA),  UK. 

The  Cases  1-5  as  classified  and  priortized  by  TTCP  workgroup  were  selected  for  the  present  study.  The 
No.  1  priority  case  which  is  the  Case  3  was  selected  for  presenting  the  results  of  grid  study.  A  standard 
grid  was  generated  in  the  Wright  Laboratory,  USA,  the  purpose  of  which  was  to  compare  results  from 
different  CFD  codes  of  all  participants  on  one  common  grid.  For  the  present  study,  the  Case  3  was  selected 
to  run  on  the  standard  grid  using  both  the  Central  and  Upwind  difference  versions.  All  the  Cases  1-5  were 
computed  for  the  most  optimum  grid  using  the  Central  difference  version  only.  Laminar  and  turbulent 
computations  were  compared  to  each  other  to  delineate  the  effects  of  the  turbulence  model  for  all  the  cases 
selected. 

Results  are  presented  for  surface  pressure  at  several  axial  stations  namely  at  X/D  of  2.4,  3.5,  4.5,  5.5, 
6.5,  7.5,  8.5  and  11.5.  Pitot  pressure  comparisons  in  the  flowfield  are  shown  depending  on  the  availability 
of  experimental  data  at  X/D  locations  of  5.5,  8.5,  and  11.5.  The  axial  and  radial  components  of  force 
and  moment  data  predicted  by  the  computational  code  for  the  standard  grid  of  Case  3  was  compared  to 
the  data  obtained  from  experiments.  Convergence  histories  for  a  laminar  and  turbulent  computation  as 
representative  of  computations  performed  to  reach  steady  state  in  the  present  study,  are  shown  in  this 
report. 

The  important  results  of  the  study  are  summarized  as  follows.  The  surface  pressure  and  pitot  pressure 
comparisons  show  very  little  variation  between  the  central  and  upwind  difference  versions  for  both  laminar 
and  turbulent  computations.  The  surface  pressure  comparisons  using  the  turbulence  model  have  matched 
the  experiment  better  than  those  with  laminar  computations.  At  an  X/D  location  between  6.5  and  8.5  at 
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(j>  between  120®  to  160®  the  surface  pressure  predicted  by  the  turbulent  computation  matches  the  laminar 
computation  and  both  do  not  match  the  experiment.  From  the  pitot  pressure  comparisons  and  the  matching 
of  the  vortex  centerlines  with  experimental  data  we  may  make  the  following  observations.  The  excessive 
dissipation  in  the  turbulence  model  in  the  capturing  of  the  vortices  is  evident  from  the  comparisons  of 
the  shape  of  the  vortices;  for  cases  1,  2,  3,  4  and  5  the  laminar  prediction  of  the  vortex  core  centerline 
matches  the  experiment  better  than  the  turbulent  prediction  at  select  X/D  locations.  The  pitot  pressure 
comparisons  show  that  the  shock  waves  outside  of  the  vortex  is  captured  in  the  turbulent  computations 
and  not  the  laminar  ones.  Below  the  primary  vortex  there  are  two  or  three  vortical  structures  for  the 
laminar  computations  and  none  for  the  turbulent  computations. 

One  of  the  main  recommendations  of  this  study  is  the  development  of  better  turbulence  models  to 
simulate  vortical  flow  features.  The  inclusion  of  turbulence  model  has  been  shown  to  be  important  in 
the  accurate  prediction  of  surface  pressure.  The  k-e  turbulence  model  has  shown  to  be  very  dissipative 
in  regions  of  vortices.  This  suggests  that  modifications  of  the  model  which  distinguishes  the  regimes  near 
wall  and  vortex.  The  other  important  area  of  development  should  be  in  the  methodology  to  improve  grid 
resolution  locally  such  as  embedded  and  adaptive  grids  to  better  capture  the  vortical  structures.  This  will 
help  to  predict  more  accurately  the  strength  and  position  of  the  vortical  flows. 
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Cased  Mach  =  2.5  AoA  =  14deg.  X/D  =  11.5 


Figure  3:  Grid  Study:  Comparison  of  Pitot  Pressure 
Figure  1:  Experimental  Data  Geometry  along  Vortex  Centerline  for  Laminar 


Grid  Used  for  Case  1:  Mach  =  1.45  AoA  =  14  deg. 


Case  3  Mach  =  2.5  AoA  =  1 4  deg.  X/D  =  1 1 .5 


Figure  2:  Computational  Grid 


Figure  4:  Grid  Study:  Comparison  of  Pitot  Pressure 
along  Vortex  Centerline  for  Turbulent 


Figure  5:  Comparison  of  Surface  Pressure  on  Stan-  Figure  8:  Comparison  of  Surface  Pressure  on  Stan¬ 
dard  Grid  at  X/D  =  2.4  for  Laminar  dard  Grid  at  X/D  =  3.5  for  Turbulent 


Figure  6:  Comparison  of  Surface  Pressure  on  Stan-  Figure  9:  Comparison  of  Surface  Pressure  on  Stan¬ 
dard  Grid  at  X/D  =  2.4  for  Turbulent  dard  Grid  X/D  =  4,5  for  Laminar 


Figure  7:  Comparison  of  Surface  Pressure  on  Stan-  Figure  10:  Comparison  of  Surface  Pressure  on  Stan¬ 
dard  Grid  at  X/D  =  3.5  for  Laminar  dard  Grid  X/D  =  4.5  for  Turbulent 
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Figure  11:  Comparison  of  Surface  Pressure  on  Stan-  Figure  14:  Comparison  of  Surface  Pressure  on  Stan¬ 
dard  Grid  at  X/D  =  5.5  for  Laminar  dard  Grid  at  X/D  =  6.5  for  Turbulent 


Figure  12:  Comparison  of  Surface  Pressure  on  Stan¬ 
dard  Grid  at  X/D  =  5.5  for  Turbulent 


Figure  15:  Comparison  of  Surface  Pressure  on  Stan¬ 
dard  Grid  at  X/D  =  7.5  for  Laminar 


Figure  13:  Comparison  of  Surface  Pressure  on  Stan-  Figure  16:  Comparison  of  Surface  Pressure  on  Stan¬ 
dard  Grid  at  X/D  =  6.5  for  Laminar  dard  Grid  at  X/D  =  7.5  for  Turbulent 
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Figure  17:  Comparison  of  Surface  Pressure  on  Stan¬ 
dard  Grid  at  X/D  =  8.5  for  Laminar 


Figure  20:  Comparison  of  Surface  Pressure  on  Stan¬ 
dard  Grid  at  X/D  =  11.5  for  Turbulent 


Figure  18:  Comparison  of  Surface  Pressure  on  Stan¬ 
dard  Grid  at  X/D  =  8.5  for  Turbulent 


Figure  21:  Comparison  of  Pitot  Pressure  on  Stan¬ 
dard  Grid  at  X/D  =  5.5  using  Central  Diff.  (Lami¬ 
nar) 


Figure  19:  Comparison  of  Surface  Pressure  on  Stan¬ 
dard  Grid  at  X/D  =  11.5  for  Laminar 
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Case  3  Mach  =  2.5  AoA  =  14  deg.  X/D  =  5.5 


Experiment  Turbulent  Computation 


(Central  Diff.  Code) 


(Upwind  DHf.  Code) 


Figure  22:  Comparison  of  Pitot  Pressure  on  stan-  Figure  24:  Comparison  of  Pitot  Pressure  on  Stan- 
dard  grid  at  X/D  =  5.5  using  Central  Diff.  (Turbu-  dard  Grid  at  X/D  =  5.5  using  Upwind  Diff.  (Lami- 
lent)  nar) 


Cases  Mach  =  2.5  AoA  =14 deg.  X/D  =  5.5 


Experiment  Turbulent  Computation 
(Upwind  Diff.  Code) 


Figure  23:  Comparison  of  Pitot  Pressure  along  Vor-  Figure  25:  Comparison  of  Pitot  Pressure  on  Stan- 
tex  Centerline  using  Standard  Grid  at  X/D  =  5.5  dard  Grid  at  X/D  =  5.5  using  Upwind  Diff.  (Tur- 
using  Central  Diff.  bulent) 
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Experiment  Turbulent  Computation 
(Central  Diff.  Code) 


Figure  26:  Comparison  of  Pitot  Pressure  along  Vor-  Figure  28:  Comparison  of  Pitot  Pressure  on  Stan- 
tex  Centerline  using  Standard  Grid  at  X/D  =  5.5  dard  Grid  at  X/D  =  11.5  using  Central  Diff.  (Tur- 
using  Upwind  Diff.  bulent) 


Experiment  ^amiiw  C^pu^lon 
(Central  Diff.  Code) 


Figure  27:  Comparison  of  Pitot  Pressure  on  Stan-  Figure  29:  Comparison  of  Pitot  Pressure  along  Vor- 
dard  Grid  at  X/D  =  11.5  using  Central  Diff.  (Lam-  tex  Centerline  using  Standard  Grid  at  X/D  =  11.5 
inar)  using  Central  Diff. 
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Case  3  Mach  =  2.5  AoA=  14  deg.  X/D  =  1 1 .5 


Figure  30:  Comparison  of  Pitot  Pressure  on  Stan-  gg.  Comparison  of  Pitot  Pressure  along  Vor- 

dard  Grid  at  X/D  =  11.5  using  Upwind  Diff.  (Lam-  Centerline  using  Standard  Grid  at  X/D  =  11.5 

using  Upwind  Diff. 


Cases  Mach  =  2.5  AoA:=14deg.  X/D  =  11.5 


Experiment  T  urbulent  Computation 
(Upwind  Diff.  Code) 


Figure  31:  Comparison  of  Pitot  Pressure  on  Stan¬ 
dard  Grid  at  X/D  =11.5  using  Upwind  Diff.  (Tur¬ 
bulent) 


Experiment  Laminar  Computation 

(Central  Diff.  Code) 


Figure  33:  Comparison  of  Pitot  Pressure  for  Case  1 
at  X/D  =  8.5  (Laminar) 
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Experiment  Turbulent  Computation 

(Central  Diff.  Code)  Experiment  Turbulent  Computation 

(Central  Diff.  Code) 


Figure  34:  Comparison  of  Pitot  Pressure  for  Case  1  Figure  37:  Comparison  of  Pitot  Pressure  for  Case  1 
at  X/D  =  8.5  (Turbulent)  at  X/D  =  11.5  (Turbulent) 


Case1  Mach  =  1.45  AoA=:14deg.  X/D=;8.5  Case  1  Mach=1.45  AoA  =  14deg.  XA:)  =  11.5 


Figure*35:  Comparison  of  Pitot  Pressure  along  Vor¬ 
tex  Centerline  for  Case  1  at  X/D  =  8.5 


Experiment  Laminar  Computation 


(Central  Diff.  Code) 

Figure  36:  Comparison  of  Pitot  Pressure  for  Case  1 
at  X/D  =  11.5  (Laminar) 


Figure  38:  Comparison  of  Pitot  Pressure  along  Vor¬ 
tex  Centerline  for  Case  1  at  X/D  =11.5 

Case  2  Mach  =  1.8  AoA=14deg.  X/D  =  5.5 


Experiment  Laminar  Computation 
(Central  Diff.  Code) 

Figure  39:  Comparison  of  Pitot  Pressure  for  Case  2 
at  X/D  =  5.5  (Laminar) 
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Case 2  Machs  1.8  AoA=14deg.  X/Ds5.5 


Experiment  Turbulent  Computation 
(Central  Diff.  Code) 


Case2  Machsl.8  AoAs14deg.  X/D  =  8.5 


Distance 


Figure  40:  Comparison  of  Pitot  Pressure  for  Case  2  Figure  43:  Comparison  of  Pitot  Pressure  along  Vor- 
at  X/D  =  5.5  (Turbulent)  tex  Centerline  for  Case  2  at  X/D  =  8.5 


Case 2  Machs  1.8  AoAsi4deg.  X/D s 8.5 


Experiment  Laminar  Computation 
(Central  Diff.  Code) 

Figure  41:  Comparison  of  Pitot  Pressure  for  Case  2 
at  X/D  =  8.5  (Laminar) 

Case 2  Machs  1.8  AoAsl4deg.  X/D  s  8.5 


Experiment  Turbulent  Computation 
(Central  Diff.  Code) 


Experiment  Laminar  Computation 
(Central  Diff.  Code) 

Figure  44:  Comparison  of  Pitot  Pressure  for  Case  2 
at  X/D  =11.5  (Laminar) 

Case2  Machs1.8  AoAs14deg.  X/D  =  11.5 


Experiment  Turbulent  Computation 
(Central  Diff.  Code) 


Figure  42:  Comparison  of  Pitot  Pressure  for  Case  2  Figure  45:  Comparison  of  Pitot  Pressure  for  Case  2 
at  X/D  =  8.5  (Turbulent)  at  X/D  =  11.5  (Turbulent) 
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Case2  Mach  =  1,8  AoA=14deg.  X/D  =  11^ 


0.0  OJt  0.4  0.0  0.0  1.0 


Distance 

Figure  46:  Comparison  of  Pitot  Pressure  along  Vor¬ 
tex  Centerline  for  Case  2  at  X/D  =  11.5 

Case  3  Mach  =  2.5  AoA  =14  deg.  X/D  =  5.5 


Experiment  Laminar  Computation 
(Central  Diff.  Code) 

Figure  47:  Comparison  of  Pitot  Pressure  for  Case  3 
at  X/D  =  5,5  (Laminar) 

Case  3  Mach  =  2.5  AoA=  1 4  deg.  X/D  =  5.5 


Experiment  Turbulent  Computation 
(Central  Diff.  Code) 

Figure  48:  Comparison  of  Pitot  Pressure  for  Case  3 
at  X/D  =  5.5  (Turbulent) 


Case3  Mach  =  2.5  AoA=14deg.  X/D  =  11.5 


Experiment  Laminar  Computation 
(Central  Diff.  Code) 


Figure  49:  Comparison  of  Pitot  Pressure  for  Case  3 
at  X/D  =11.5  (Laminar) 

Case3  Mach  =  2.5  AoA  =  14deg.  X/D  =  11.5 


Experiment  Turbulent  Computation 
(Central  Diff.  Code) 


Figure  50:  Comparison  of  Pitot  Pressure  for  Case  3 
at  X/D  =  11.5  (Turbulent) 


Case  3  Mach  =  2.5  AoA  s  14  deg.  X/D  =  11 .5 


Distance 


Figure  51:  Comparison  of  Pitot  Pressure  along  Vor¬ 
tex  Centerline  for  Case  3  at  X/D  =  11.5 
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Case  4  Mach  =  3.5  AoAsSdeg.  X/D  =  11.5 


Figure  52:  Comparison  of  Pitot  Pressure  for  Case  4 
at  X/D  =  11.5  (Laminar) 


Cases  Mach  =  3.5  AoA=14deg,  X/D  =  5.5 


E^qTeriment  Laminar  Computation 
(Central  Diff.  Code) 

Figure  55:  Comparison  of  Pitot  Pressure  for  Case  5 
at  X/D  =  5.5  (Laminar) 

Cases  Mach  =  3.5  AoA  =  14deg.  X/D  =  5.5 


Experiment  Turbulent  Computation 
(Centra!  Diff.  Code) 


Figure  53:  Comparison  of  Pitot  Pressure  for  Case  4  Figure  56:  Comparison  of  Pitot  Pressure  for  Case  5 
at  X/D  =  11.5  (Turbulent)  at  X/D  =  5.5  (Turbulent) 


Case4  Mach  =  3.5  AoA  =  8deg.  X/D  =  11.5 


Figure  54:  Comparison  of  Pitot  Pressure  along  Vor¬ 
tex  Centerline  for  Case  4  at  X/D  =  11.5 


Cases  Mach  =  3.5  AoA  =  14deg.  X/D  =  11.5 


Experiment  Laminar  Computation 


(Central  Diff.  Code) 

Figure  57:  Comparison  of  Pitot  Pressure  for  Case  5 
at  X/D  =11.5  (Laminar) 
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Figure  58:  Comparison  of  Pitot  Pressure  for  Case  5  Figure  61:  Comparison  of  Surface  Pressure  for  Case 
at  X/D  -  11.5  (Turbulent)  1  at  X/D  =  3.5 


Cases  Mach 3.5  AoA=14deg.  X/D  =  11.5 


Case  1  Mach  =  1 .45  AoA  =  14  deg.  X/D  =  4.5 


Figure  59:  Comparison  of  Pitot  Pressure  along  Vor-  Figure  62:  Comparison  of  Surface  Pressure  for  Case 
tex  Centerline  for  Case  5  at  X/D  =  11.5  1  at  X/D  =  4.5 


Case  1  Mach  =  1 .45  AoA  =14  deg.  X/D  =  2.4 


Case  1  Mach  =  1 .45  AoA  =  14  deg.  X/D  =  5.5 


Figure  60:  Comparison  of  Surface  Pressure  for  Case  Figure  63:  Comparison  of  Surface  Pressure  for  Case 
1  at  X/D  =  2.4  1  at  X/D  =  5.5 


Case1  Mach  =  1.45  AoA=  14  deg.  X/D  =  6.5  Case  1  Mach  =  1.45  AoA  =  14 deg.  X/D  =  115 


Figure  64:  Comparison  of  Surface  Pressure  for  Case  Figure  67:  Comparison  of  Surface  Pressure  for  Case 
1  at  X/D  =  6.5  1  at  X/D  =  11.5 


Figure  65:  Comparison  of  Surface  Pressure  for  Case  Figure  68:  Comparison  of  Surface  Pressure  for  Case 
1  at  X/D  =  7.5  2  at  X/D  =  2.4 


Figure  66:  Comparison  of  Surface  Pressure  for  Case  Figure  69:  Comparison  of  Surface  Pressure  for  Case 
1  at  X/D  =  3.5  2  at  X/D  =  3.5 
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Figure  70:  Comparison  of  Surface  Pressure  for  Case  Figure  73:  Comparison  of  Surface  Pressure  for  Case 
2  at  X/D  =  4.5  2  at  X/D  =  7.5 


Phi 


Case2  Mach=:1.8  AoA  =  14deg.  X/D-8.5 


Figure  71:  Comparison  of  Surface  Pressure  for  Case  Figure  74:  Comparison  of  Surface  Pressure  for  Case 
2  at  X/D  =  5.5  2  at  X/D  =  8.5 


Case2  Mach=:1.8  AoA  =  14deg.  X/D  =  11.5 


Figure  72:  Comparison  of  Surface  Pressure  for  Case  Figure  75:  Comparison  of  Surface  Pressure  for  Case 
2  at  X/D  =  6.5  2  at  X/D  =  11.5 
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Figure  76:  Comparison  of  Surface  Pressure  for  Case  Figure  79:  Comparison  of  Surface  Pressure  for  Case 
3  at  X/D  =  2.4  3  at  X/D  =  5.5 


Figure  77:  Comparison  of  Surface  Pressure  for  Case  Figure  80:  Comparison  of  Surface  Pressure  for  Case 
3  at  X/D  =  3.5  3  at  X/D  =  6.5 


Figure  78:  Comparison  of  Surface  Pressure  for  Case  Figure  81;  Comparison  of  Surface  Pressure  for  Case 
3  at  X/D  =  4.5  3  at  X/D  =  7.5 


35 


Case  3  Mach  =:  2.5  AoA  =14  deg.  XTD  =  8.5 


Figure  82:  Comparison  of  Surface  Pressure  for  Case  Figure  85:  Comparison  of  Surface  Pressure  for  Case 
3  at  X/D  =  8.5  4  at  X/D  =  3.5 


Cases  Mach  =  2.5  AoA  =  14deg.  X/D  =  11.5 


Figure  83:  Comparison  of  Surface  Pressure  for  Case  Figure  86:  Comparison  of  Surface  Pressure  for  Case 
3  at  X/D  =  11.5  4  at  X/D  =  4.5 


Phi 


Case  4  Mach  =  3.5  AoA  =  8  deg.  X/D  =  5.5 


Phi 


Figure  84:  Comparison  of  Surface  Pressure  for  Case  Figure  87:  Comparison  of  Surface  Pressure  for  Case 
4  at  X/D  =  2.4  4  at  X/D  =  5.5 
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Case4  Mach=:3^  AoA  =  8deg.  XA}-11.5 


Phi 


Figure  88:  Comparison  of  Surface  Pressure  for  Case  Figure  91:  Comparison  of  Surface  Pressure  for  Case 
4  at  X/D  =  6.5  4  at  X/D  =  11.5 


Figure  89:  Comparison  of  Surface  Pressure  for  Case  Figure  92:  Comparison  of  Surface  Pressure  for  Case 
4  at  X/D  =  7.5  5  at  X/D  =  2.4 


Case  4  Mach  =  3.5  AoA  =  8deg.  X/D  =  8.5 


Phi 


Figure  90:  Comparison  of  Surface  Pressure  for  Case  Figure  93:  Comparison  of  Surface  Pressure  for  Case 
4  at  X/D  =  8.5  5  at  X/D  =  3.5 
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Case  5  Mach  =  1 .8  AoA  =14  deg.  X/D  =  4.5  Case  5  Mach  =  3.5  AoA  =14  deg.  X/D  =  7.5 


Figure  94:  Comparison  of  Surface  Pressure  for  Case  Figure  97:  Comparison  of  Surface  Pressure  for  Case 
5  at  X/D  =  4.5  5  at  X/D  =  7.5 


Case  5  Mach  =  3.5  AoA  =  14  deq.  X/D  =  5.5  Case  5  Mach  =  3.5  AoA  =  14  deg.  X/D  =  8.5 


Figure  95:  Comparison  of  Surface  Pressure  for  Case  Figure  98:  Comparison  of  Surface  Pressure  for  Case 
5  at  X/D  =  5.5  5  at  X/D  =  8.5 


Cases  Mach  =  3.5  AoA-14deg.  X/D  =  6.5  CaseS  Mach-3.5  AoA  =  14deg.  X/D  =  11.5 


Figure  96:  Comparison  of  Surface  Pressure  for  Case  Figure  99:  Comparison  of  Surface  Pressure  for  Case 
5  at  X/D  =  5.5  5  at  X/D  =  11.5 
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Figure  100:  Convergence  History  for  Laminar  Com¬ 
putation 


Figure  101:  Convergence  History  for  Turbulent 
Computation 
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